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Abstract 

We present a model for a thermal Brownian motor based on Feynman's famous 
ratchet and pawl device. Its main feature is that the ratchet and the pawl are in 
different thermal baths and connected by an harmonic spring. We simulate its dy- 
namics, explore its main features and also derive an approximate analytical solution 
for the mean velocity as a function of the external torque applied and the temper- 
atures of the baths. Such theoretical predictions and the results from numerical 
simulations agree within the ranges of the approximations performed. 
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1 Introduction 



The engines which aim to get useful work by rectifying thermal fluctuations 
are called Brownian motors (BM). During the last years a lot of effort has been 
invested to study the underlying mechanism of such engines which has been 
called as the ratchet effect [1] . This is a mechanism which consists in breaking 
the spatial and temporal inversion symmetry of the system so that directed 
transport emerges, being thermal fluctuations the very relevant input. In fact, 
the paradigmatic device of such speculations is Feynman's famous ratchet and 
pawl machine [2]. 

In 1963 R. P. Feynman introduced [2] a microscopic device (the ratchet and 
pawl machine) that can operate between two thermal baths extracting some 
mechanical work. The hotter bath contains an axle with vanes in it. The 
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Figure 1. Feynman's original ratchet and pawl machine. 



bombardments of gas molecules on the vane make the axle rotate with random 
symmetric fluctuations. At the other end side of the axle, as shown in Fig. (1), 
there is a second box with an asymmetric toothed wheel which in principle can 
turn only one way due to a coupling with a pawl (the stopping mechanism). 

At first glance one could think that it seems quite likely that the wheel will 
spin round one way and lift a weight even when both gases are at the same 
temperature and thus violating the Second Law. However, a closer look at the 
pawl reveals that it bounces and so the wheel will rotate randomly in any 
direction, doing a lot of jiggling but with no net turning. Thus, the machine 
cannot extract work from two baths at the same temperature. When the tem- 
perature of the vanes is higher than the temperature of the wheel, Feynman 
concludes that some work is performed with Carnot's efficiency when the ma- 
chine is lifting the weight very slowly. This is indeed a very optimistic result 
which has been revised, many years later, in Refs. [3,4]. Such particular device 
has been analyzed and it has been suggested that this engine is very far from 
Carnot efficiency [3] due to the fact that this type of device has strong heat 
losses. A more refined analysis and modeling of Feynman's motor has been 
presented in several references [4,5,6], with the conclusion that its efficiency 
is very poor. The main problem comes from the mechanical coupling of the 
pawl mechanism, which is not very efficient. 

Moreover, the main idea underlined in the the ratchet and pawl mechanism 
can be implemented in different ways. In this work we propose a model through 
an equation for a dynamical classical variable (position or angle) moving in a 
periodic and asymmetric potential (a ratchet potential) in a bath at tempera- 
ture Ti coupled with another degree of freedom immersed in a different bath at 
temperature T2. The coupling mechanism is a harmonic spring. In this case the 
two thermal baths are clearly separated and the performance of this engine, as 
a function of the different parameters, can be studied. In the following section, 
we introduce the model. Once the model is explained and justified, we present 
in Section 3 results by numerical simulations of the stochastic equations of 
motion. In Section 4 an analytical approach is presented and in Section 5 we 
show the agreement of the theoretical expressions found with the data from 
numerical simulations. Finally we end with some conclusions and comments. 
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Figure 2. Ratchet, Pawl and Spring Brownian Motor. 
2 The Ratchet, Pawl and Spring Motor Brownian Motor (RPSBM) 

First of all, let's describe our proposal for a Brownian motor as it is shown in 
Fig. (2). It consists of two boxes at different constant temperatures. The left 
box contains a ratchet and a pawl that act like a mechanical rectifier device. 
From this box an axel-wheel device is connected allowing it to lift a hanging 
object in order to study how much work the device can perform. The second 
box is at a higher temperature than the first one and has a little windmill that 
is used to pick up energy from the thermal bath (through the collisions of the 
particles of the bath with the vanes) and to transfer it to the ratchet and pawl 
system through the spring. The main simplification of our model is that the 
dented wheel and the pawl mechanisms are substituted by a ratchet potential 
and a harmonic spring. Then the Langevin equations in the over-damped limit 
are straightforward to write, 

rlR 

)^2^ = m-e2) + Ut), (2) 



with thermal noise satisfying the fluctuation-dissipation relation, 

iUtMt')) = 2kBTiXAjS{t - t'). (3) 



Note that we are describing the evolution in time of the angular positions 6i 
and 02 of the ratchet and the windmill respectively. We see in equation (1) 
a force term that comes from a periodic and asymmetric potential Vr that 
models the shape of the sawtooth wheel. There's also an external torque r 
that can be different than zero to account for the work done. Two independent 
white noises and ^2 represent thermal fluctuations and, flnally, there's an 
interaction between the two baths through a harmonic spring of constant k. 
A hnear coupling between both degrees of freedom also appears in Ref . [7] . 
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Figure 3. Shape of the ratchet potential Vr ioi d = 4. 
The ratchet potential Vr(^) is given by, 



2.23 



[sin(ci^) + 0.275 sin(2(i^) + 0.0533 sin(3ci^)], 



(4) 



where Vq controls the height of the potential, d is the number of teeth per turn, 
and the asymmetry of the potential is controlled by changing the numerical 
coefficients that multiply the sinus functions. In Fig. (3) we see the explicit 
form of the ratchet potential Vji{9) used in our study. 

This problem has a set of parameters which can be simplified if a new time 
dimensionless scale is defined, t = y^s. Then the equations (1) and (2) are 
transformed and all the relevant parameters: Ti = -^p^, T2 = ^y'^^ , ^ = it^j 
A = ^ and f = y^, are now dimensionless. In this way we can concentrate on 
the main parameters that control the dynamics of the system. For instance, 
we see that Vq controls the energy scale, with respect to the thermal energy 
of the baths. 



3 Numerical results and optimal regime of the motor 

In this section we present the results by numerically simulating Eqs. (1) and 
(2). We have used a second order algorithm (Heun), which is a generaliza- 
tion of the Runge-Kutta algorithm for stochastic systems. Since the process 
is intrinsically nondeterministic, it is convenient to do statistics in order to 
minimize fluctuations in the output data. For that reason, we have averaged 
every step of the integration over 500 different particles. In figure (4) we show 
an example of a single trajectory in time and the mean (^i(s)). 

We must first explore the role played by every parameter of the model in order 
to get in scale and tune them cleverly to obtain the most efficient working 
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Figure 4. Typical evolution of the angle position 6i as a function of the dimensionless 
time s directly obtained from the simulations. The parameters are Vq = SksTi, 
k = lOOksTi, d=l2 and T2 = 2Ti. Plot of a single trajectory (for r = 0) in which 
one can see the discrete jumps of the dynamics from one valley to another of the 
saw-tooth potential. Inset: Average of ^i(s) over many particles for r = O.OlkBTi. 

regime. The height of the potential Vq is a very relevant parameter since it 
determines the energy scale of the system. We are interested in an energy 
barrier higher than the energies of the two thermal baths but not much more. 
Otherwise, the engine wouldn't work because thermal fluctuations would not 
be able to cause a jump over such a high barrier. There is an optimal value 
that is found to be around Vq = 3KbTi. For smaller values, the height is 
too small and thus fluctuations can jump it over to both sides, giving a zero 
velocity on average. For large Vq, hardly ever is the system able to perform a 
jump to the next valley. This is what is expected from intuition and also what 
is found from numerical simulations. 

Regarding temperatures, we must choose T2 higher than Ti but not too much 
because large randomness [8] blurs and destroys directional motion. In a fol- 
lowing section we will analyze in detail the behavior of the mean speed as a 
function of the ratio of temperatures. In fact, there is an optimal ratio which 
is around T2 — 4Ti. 

We have also checked that the velocity falls to zero for two limits of the 
coupling constant: k and k ^ (X). This is in perfect accordance with the 
following thermodynamical arguments. For k going to zero, the systems are 
un-coupled and therefore we have a single bath at a constant temperature. 
In this situation the wheel will have no mean velocity. For the opposite limit, 
this is, large values of the coupling constant, we use the argument discussed 
in Ref. [3]. A very high k is equivalent to having both systems joined by a 
rigid axel vane. Then, the system is in contact with two thermal baths at 
different temperatures. One can prove that the system feels a single averaged 
temperature and thus it will not run at all on average either. 
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There seems to be an optimal number of teeth d as well. For low values of d 
the wheel isn't a proper sawtooth and the asymmetry of the potential is not 
felt. For d around d = 24. we find that the wheel spins quicker, while for d 
greater the velocity decays again. Notice that the force term is proportional 
to the number of teeth. This means than the higher d is, the steepest the 
potential becomes. 

Using three other different models of the ratchet potential, which are built 
similarly but modifying the coefficients multiplying the sinus functions, we see 
that the more asymmetrical the model is, the faster the system runs. Obviously 
for a symmetric potential no rectification of the fiuctuations is possible for any 
value of the rest of parameters. 

When examining equations (1) and (2) in dimensionless variables, one finds 
that the noises, in dimensionless variables too, should be written as 

{Us)Us')) = 2T\S{s-s'), (5) 

{i2{sMs')) = 2XT2S{s-s'), (6) 

{ii{s)i2{s'))^0. (7) 



Notice that the effective intensity of the second bath is the thermal energy of 
the bath in terms of Vq multiplied by the fraction of the friction coefficients. 

Finally, when setting Ai = A2, = 2Ti, Vq = 3KbTi, t = 0, k ^ SOksTi and 
d = 24, we are in an optimal regime in which the an angular velocity is nearly 
the greatest possible: v = (^1) ~ 0.021s~^. 

Apart from the quantitative numerical value itself, we can make two main 
conclusions. The first one is that the engine does work. We have made some 
simulations in which we use a torque r different than zero. For instance, for 
a small torque r = OAksTi we see that the systems still runs and lifts the 
external weight, thus performing useful work. The second conclusion is that 
this motor is very inefficient. Energetics are not worth to be calculated in detail 
because just by comparison to the model in [5] we can see that the maximum 
velocity achieved by the RPSBM is at least five times smaller than the speeds 
found in Ref. [5] for the SBM. Therefore the efficiency of the present model is 
even much smaller. 



4 Analytical study of the RPSBM 



Our purpose now is to make an analytical study of our motor. Any kind 
of exact calculation in such systems seems impossible and, therefore, some 
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approximations have to be assumed. 



Consider the equations that define our Brownian motor written in the form 

^■l = /(^l)-^(^l-^2)-T + 6W, (8) 

e2 = k{ei-e2) + ut), (9) 

where f{Oi) is the force exerted by the ratchet potential: f{9i) = —V^{9i). 
Note that we have set the friction coefficients to one for simphcity without 
losing any generality in the calculus. Let us define now the changes of variables 

n ' n ■ (10) 



The relevant variable x describes the evolution of the center of mass and 
the "irrelevant" variable y describes the relative motion of the two-particle 
system. The equations of motion in these new variables are 

i = /(£±»)-I + „M, (11) 

. ^ /(£ + ») _^_2ky + ,,((), (12) 



where a redefinition of the noises has been introduced as 

. iiW±M). . iiW^. (13) 



Let's make our first approximation. Using the fact that y is very small, we can 
make a Taylor expansion of /(a; + y) up to first order in y, obtaining the pair 
of equations 

,^/M±i^M^ + (14) 
i = M^i±Mm^ + ^(,). (15) 



Since numerical simulations indicate that the variable y has a faster dynamics 
we will eliminate it adiabatically. This means we set ^ = 0. Then Eq. (15) 
reduces to 
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One can see that the term 4k is much bigger than f'{x) when k ~ lOOksTi 
and d is small (d ~ 4). Then, within this range of parameters we can keep only 
4/c in the denominator. Substituting now the last expression for y in equation 
(14) we end up with a Langevin equation with a new force H{x) and two 
multiplicative noises gi{x) and g2ix), 



X = H{x) + gi{x)Ut) + 92{x)Ut). (17) 



where, 



H{x) = \{f{x)-T){l+^-^l 



(18) 



The Fokker-Planck equation associated to Eq. (17) is 

dtP{x,t) ^ -d,J{x,t), (20) 

where 

J{x,t) = H{x)P{x,t) - kBT,[g,{x)d^gi{x)P{x,t)] 

-kBT2[g2ix)d,g2{x)P{x,t)]. (21) 

After some manipulations with partial derivatives, the probability current 
J{x, t) can be rewritten as 

J{x,t) = H{x)P{x,t) - [geff{x)d,geffix)P{x,t)], (22) 

with 



geffix) = ^JkBT^gf{x)+kBT2gi{x). (23) 



Since we arc interested in the steady state, we have to solve now equation (22) 
for a constant probability current J. The first step is to reduce this equation 
to a Bernoulli form which can be formally integrated. By imposing periodic 
boundary conditions, P{x) = P{x + L), where L = ^, we find that 

Po(l - e^(^)) = J / dx --, (24) 
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where Pq is a constant that can be found from the normahzation condition 
J^P{x)dx - 1, and p is a relevant function whose expression is 



The mean velocity v = (x) is just the integral of the probability current over 
the spatial period L. Since J is constant, we have from Eq. (24) that the 
solution for v is [9,10,11,12,13] 

v^J^fil-ef"), (26) 



where Af is a, constant that is found by using the normalization condition for 
the probability P{x). We will see that this constant depends very smoothly on 
the control parameters of our model. Notice that the mean velocity v ~ (^i) 
because (6'i) ~ {62). 

At this point, what is left to do is to find the /3 integral as a function of the 
parameters we are interested in. From now on we will assume implicitly that 
the energy is in ksTi units. To simplify the calculation of the integral, we 
make an expansion of the denominator in equation (25) in powers of |. Since 
the value of k that we will consider is the one that makes the motor run faster 
{k ~ 100), one can safely suppose that the terms of the order (|)^ and so on 
will not notably contribute to the integral. Remember that d and Vq are kept 
small. The relevant quantity /? can then be written as 

^ / ill + 1 + 2(1 -t)^l- '^'^ 



Let's note that the terms ~f{x) + r in the numerator are very small when 
integrated. Therefore we can neglect the much smaller correction f'[x)/4:k 
of the numerator. However, one cannot do this approximation for the same 
term in the denominator because it is the responsible from the net motion of 
the Brownian motor. Such term contains two essential physical features. The 
first one is that it accounts for the multiplicative noise and thus, without it, 
thermal sources are unable to make the motor move. The second one is that 
it avoids the violation of the Second Law, i.e. when T2 = Ti it cancels and no 
average velocity is predicted. Then, the final and simplest expression for /3 is 

-fix) + T 



/5 = 2 / ^ ' ^ dx. (28) 

' 1 + ^ + 2(1-1)^ 



9 



0,0005 
0,0004 
0,0003 
0,0002 
0,0001 


0,001 0,002 0,003 0,004 0,005 
T 

Figure 5. Comparison between numerical simulations (dots) and theoretical predic- 
tion (line) from (33) of the dependence of the mean velocity v as a function of the 
external torque applied r. The parameters chosen are k = 100, d = A, Vq = 3 and 
T2 = 2Ti. 




5 Simulations versus theory 



Now that we have derived an analytical expression for the mean velocity of 
the motor, let's compare its predictions to the numerical simulations when 
exploring two relevant parameters: the external torque r and the temperature 
difference ratio T2/T1. 

The first comparison is plotted in figure (5), in which we show the mean ve- 
locity as a function of the external torque. It is very remarkable that the stall 
torque (the torque at which the motor stops) is perfectly predicted by the ex- 
pression found. What is more, the functional behavior (the linear dependence 
with negative slope) is clearly reproduced. To determine the scale we adjust 
J\f so that the analytical prediction fits the simulations at r = 0, at the same 
time that we confirm that such constant J\f does not depend very much on r. 

Since (3 is found to be very small, one can write the following expression for 
the mean velocity, 

v^Af{l-e>^)^Af{-p)^a-rb, (29) 



where 
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Figure 6. Mean velocity v as a function of the ratio of temperatures T2/T1. The 
analytical prediction (line) fits the numerical results (dots) quite well allowing for 
the fact that there's some uncertainty in the data obtained from the simulations. 
In this case r = and the other control parameters are unchanged. 

) 1 

b^2Af —^rrrdx. (31) 

J 1 + I2 +2(1 - 2i)£M ^ ^ 



This last expression for the parameter b can be easily simplified by expanding 
its denominator and stopping at the first term, finding 

bc^— ^. 32 



The term a is a little more complicated but it can be computed numerically. 
Finally we obtain 

V ~ 0.00049 - 0.097r, (33) 



in units of s ^. This result is plotted in Fig. (5). 

Our second result is shown in figure (6), in which the prediction for the mean 
velocity is v = a, as one finds when setting r = in Eq. (29). Despite data 
precision difficulties in the simulations, we can say that theory and simulations 
agree quantitatively and also reproduce three important qualitative facts. The 
first is that at 71 = T2 the mean velocity is zero, as we explained before. The 
second is that there is a maximum of the speed (an optimal value for the ratio 
of temperatures). Thirdly, we observe a slow decay of the velocity to zero as T2 
is increased. We underline as well that the constant A/" used now is the same 
that we found for the v{t) plot in Fig. (5). This again confirms our hypothesis 
that J\f depends very smoothly on the parameters we are exploring (r and 
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T2/T1). When expanding the denominator of equation (30), we find that the 
terms appear as a function of the ratio T2/T1 and its powers. Therefore, we 
can roughly see that this more complex dependence of v{T2/Ti) comes from 
such terms. However, one must keep in mind that it is not a good idea to try to 
find exactly which these terms are since our very first approximations killed 
terms of the same order in k. Nevertheless, we can say that our analytical 
predictions reproduce very well the numerical results. They capture the main 
qualitative behavior and also fit the quantitative data. 



6 Conclusions 

We have presented and studied a model for a thermal Brownian motor inspired 
on Feynman's ratchet and pawl [2]. After explaining Feynman's idea we have 
introduced and justified the RPSBM model and its main features have been 
explored numerically, finding the optimal regime of the motor. The results of 
the simulations are consistent with fundamental physical arguments that must 
always hold. We have performed an analytical calculation based on [5] with 
appropriate approximations to get an expression for the mean velocity in terms 
of the relevant parameters of the model. We have analyzed its dependence 
on the external torque r and the ratio of temperatures T2/T1. Such formal 
predictions fit very well the data from the numerical simulations. 

Regardless of the particular properties of these kind of heat engines, they 
are anyhow unrealistic models for molecular motors [14] since it is known 
that such biological systems are mono-thermal and convert chemical energy 
into work, without the intermediate state of burning fuel. Consequently, one 
cannot think of these models as realistic ones for biological molecular motors. 
Moreover, the mechanical coupling mechanism between both baths acts as 
a very good heat conductor even in situations of very small mean velocity. 
Therefore, the efficiency is only a small fraction of that of Carnot. What is 
more, it is much smaller than the efficiency found in similar models [5] , as one 
can see from the values of the velocity v. Such very low efficiency is a general 
feature of Brownian motors due to the fact that, in order to rectify thermal 
fluctuations, these systems must be tightly connected and then a lot of heat 
is interchanged. 

This research was supported by the Ministerio de Educacion y Ciencia (Spain) 
under project BFM2003-07850-C03-01 and the Generalitat de Catalunya (Spain) 
under grant 2005FI 00194. 
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